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The global rotational degrees of freedom in the Schrodinger equation 
for an A^-body system are completely separated from the internal ones. 
After removing the motion of center of mass, we find a complete set of 
(2£ + 1) independent base functions with the angular momentum I. These 
are homogeneous polynomials in the components of the coordinate vectors 
and the solutions of the Laplace equation, where the Euler angles do not 
appear explicitly. Any function with given angular momentum and given 
parity in the system can be expanded with respect to the base functions, 
where the coefficients are the functions of the internal variables. With the 
right choice of the base functions and the internal variables, we explicitly 
establish the equations for those functions. Only (3A^ — 6) internal variables 
are involved both in the functions and in the equations. The permutation 
symmetry of the wave functions for identical particles is discussed. 
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I. INTRODUCTION 

For a quantum A^-body system with a pair potential, the Schrodinger equation is 
invariant under spatial translation, rotation, and inversion. It is well known that, due 
to the translation symmetry of the system, the wave function can be separated into a 
product of two parts. One describes the motion of the center of mass as a free particle, 
and the other describes the motion of the system in the center-of-mass frame. It is no 
loss of generality to suppose the center of mass of the system to be at rest, so that the 
configuration is completely specified by (A^ — 1) vectors Tcj, 1 < j < — 1, which are 
usually chosen as the Jacobi coordinate vectors Rj for simplicity (see Sec. II). On 
the other hand, due to the symmetries of the global rotation and space inversion of the 
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system, the three rotational degrees of freedom should be separated completely from the 
internal ones so that only (3A^ — 6) internal variables, called the shape coordinates in 
some papers, are involved both in the functions and in the equations. This is the aim of 
this paper. 

The hydrogen atom is a typical quantum two-body system, where there is only one 
Jacobi coordinate vector, which is proportional to the relative position vector r = ri — r2. 
The Schrodinger equation for the hydrogen atom becomes a partial differential equation 
with respect to three components of r. Because of the spherical symmetry, the angular 
momentum is conserved, and the wave function can be expressed as a product of a radial 
function 0(r) and a spherical harmonic function Y^{9,Lp), 

¥^{v)=<P{r)Yi{e,v), (1) 

where the radial function (j){r) satisfies the radial equation containing only one radial 
variable. The generalization of this method to a quantum A^-body system is an important 
and fundamental problem that has been attacked by many groups. 

Wigner studied this problem using group theory Q]. The global rotation of a system 
can be described by a space rotation R = 7), rotating the center-of-mass frame 

into the body-fixed frame, where a, f3 and 7 are the Euler angles. Briefly denoting all the 
internal variables for simplicity by ^, which is invariant in the global rotation, one may 
express the wave function with a given angular momentum as \E'^(i?, ^) = \E'f^(a, /5, 7, ^). 
Let Ps be the transformation operator for a scalar function iIj{x) in the transformation 
S, Psip{x) = ip^S^^x) (see p. 105 in Ref. 0). In a rotation S G S0(3), the function 
\E'^(i?, ^) transforms as 

m'=-e 

\E'^(-R, ^) was called by Wigner the function belonging to the mth row of the repre- 
sentation D^(S0(3)). Letting R = R{0, 0, 0) be the identity element, one obtains 

m'=-l 

where ^^'(0 = ^In'^i 0; 0; depends only on the internal variables, called the generalized 
radial functions in this paper. Due to the spherical symmetry, one only needs to study the 
eigenfunctions of angular momentum with the largest eigenvalue of (m = £), which in 
this paper are simply called the wave functions with angular momentum i for simplicity. 
Their partners with the smaller eigenvalues of Lz can be calculated from them by the 
lowering operator L_. Letting S^^ = -R(a,/?, 7), one obtains (see Eq. (19.6) in f^) 

vi/^(a,A7,0 = ^^(^,0= E DU<y,p,jr4>LiO- (2) 

m=—£ 
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where the commonly used form of the ^function P is adopted. In Eq. (2) Df^{a, 
plays the role of the base function. What Wigner proved is that {2i + 1) functions 
Df^{a, P,'-yy constitute a complete set of independent base functions with the angular 
momentum £, and any wave function with the angular momentum i can be expanded with 
respect to those base functions. Due to the singularity of the Euler angles, the generalized 
radial equations satisfied by the generalized radial functions are quite difficult to derive 
based on Eq. (2). Hirschfelder and Wigner studied the problem of the generalized radial 
equations. Later, the generalized radial equations were improved by several authors 0-0. 
The equations seem quite cumbersome [0 . In the present paper we are going to rechoose 
the base functions as the homogeneous polynomials in the components of the coordinate 
vectors so that the derivation of the generalized radial equations becomes very simple. 

The generalized radial equations for a quantum three-body system have been discussed 
in more detail [p^|-[T^. Recently, by making use of the body-fixed frame, the expression 
for the kinetic energy operator was built in terms of the partial angular momentum op- 
erators and radial derivatives containing (3iV — 6) internal variables []T3|,|T1|. A coupled 
angular momentum basis was used to prediagonalize the kinetic energy operator, where 
some off-diagonal elements remain nonvanishing. Those results have been generalized to 
nonorthogonal vectors [T3,rB|. In those calculations, a function with a given angular mo- 
mentum was obtained from the partial angular momentum states using Clebsch-Gordan 
coefficients. Since the partial angular momenta are generally not conserved, one has to 
deal with, in principle, an infinite number of partial angular momentum states. This 
problem also occurs in the hyperspherical harmonic function method and its improved 
versions [0-|22|. It causes unnecessary degeneracy of the hyperspherical harmonic states 



because, as Wigner proved, only [2i + 1) base functions with angular momentum i are 
involved in the calculation. 

Eckart IQ presented another method, called the principal axis transformation, to dis- 
tinguish the global rotation and the internal motion in a classical A^-body system. From 
the (A^ — 1) Jacobi coordinate vector Rj, he defined an (A^ — 1) x (A^ — 1) real symmetric 
matrix TZji^ = Hj ■ R^, which is semipositive definite. TZ can be diagonalized by a real 
orthogonal similarity transformation X, XTZX~^ = T. 

yYl ^js^sj ■ yYl ^kt^tj =SjkTjj. (3) 

Since there are at most three orthogonal vectors in a three dimensional space, the vectors 
can be expressed as 

Y: X,,K, = h^'^ whenj = a<3 

t^i [0 when 4 < J < A^ - 1, ^ ^ 

where are three orthonormal vectors in the usual three dimensional space and = Taa- 
Three orthonormal vectors contain three Euler angles describing the global rotation of 
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the system M. Thus, from Eqs. (3) and (4), Eckart obtained 



Rj=E ^araXaj, 1<J<N-1. 



(5) 



Xaj are the first three rows of the matrix X and contain (3A^ — 9) independent param- 
eters. Therefore, VaXaj contain (3A^ — 6) internal variables, and the rotational variables 
(Euler angles) can be completely separated in the kinetic energy expression from the 
internal variables This approach has been further studied and quantized in recent 
years P,p3|-25|. The internal coordinates and their conjugate momenta were quantized to 
derive the kinetic energy expression through generalized angular momentum operators. 
However, the formal formula (5) does not give the explicit functional relation of the inter- 
nal coordinates with the components of the Jacobi coordinate vectors, so that the kinetic 
energy expression cannot be transformed directly from the usual kinetic energy term in 
the Schrodinger equation by replacement of variables. It is very difficult to obtain the 
wave function on the position vectors (or on the Jacobi coordinate vectors Rj) from a 
solution on these internal coordinates. The intermediate calculations for the kinetic en- 
ergy expression are so complicated that, as said in Ref. ||2^, the expression for a quantum 
six-body system has not been obtained probably due to a few mistakes in calculations. 

Let us return to the hydrogen atom problem. After removing the motion of the center 
of mass, the configuration space is parametrized in terms of the rectangular coordinates 
r = {x,y,z) or the spherical coordinates {r,6,ip), where r specifies the internal (radial) 
motion and {6, (p) specify the overall rotation. There is another way to separate the 
rotational degrees of freedom and obtain the same radial function and the radial equation 
as those derived from Eq. (1). One may avoid explicitly introducing the rotational angles 
9 and (p using the harmonic polynomial y^i'^) = which is a homogeneous 

polynomial of degree £ in the rectangular coordinates (x, y, z) and satisfies the Laplace 
equation as well as the eigenequation of the angular momentum. Using y^{v), Eq. (1) 
becomes 



Under the action of the Laplace operator, we have 



A^l.(r) = yl{v 

= yLi^y-'dir 



Ar" 



+ 2v[r-V(r)]-v{3^^(r)} 
+ 29, [r-V(r)] r-^r- v{3^i(r)} 



yl{Y)r~^ [r-^dlr(j){r) + ^{^ + l)r-^(j){r) + 2 (-fr^^) [9,0(r) + r" V(r)] } 
+2 [-^r-'~^<P{r) + r-'dr<p{r)] {^r~^yl(v)] 
= Yiie, if) {r-i9>0(r) - £{i + l)r-20(r)} , 
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where and hereafter dril' denotes dip/dr and so on. The results are the same. In the tradi- 
tional approach, the property that Y^{0, ip) is the eigenfunction of is used, and in the 
new approach, the property that 3^m(^) ^ homogeneous polynomial in the rectangular 
coordinates and the solution to the Laplace equation is used. It is worth emphasizing 
that the rotational angles 6 and do not appear explicitly in this approach. Since the 
differential calculus with respect to 9 and (p is not complicated, this approach is similar 
to the traditional one in a two-body system. However, it may be easier in an A^-body 
system due to the complicated calculus with respect to the Euler angles. In the present 
paper we will separate the global rotational variables in the Schrodinger equation for an 
A^-body system from the internal ones by a generalized method following the approach 
described above. In our approach, the number of base functions with the given angular 
momentum is finite, but that number in the hyperspherical harmonic function method 
and its improved versions |lT3|JT7| -p2| is infinite due to the unconserved partial angular 
momenta. We also avoid the heavy differential calculus with respect to the Euler angles 
which is sometimes necessary for expressing kinetic energy operators. 

This paper is organized as follows. In Sec. II we will briefly review the method of 
separating the motion of center of mass by the Jacobi coordinate vectors. In Sec. Ill we 
will define (3A^ — 6) internal variables from the Jacobi coordinate vectors Rj and find the 
{2^+ 1) base functions with total orbital angular momentum which are the homogeneous 
polynomials in the components of Rj and the solutions of the Laplace equation. Then we 
will prove that the base functions constitute a complete set, namely, any function with 
the angular momentum £ and the given parity in the system can be expanded with respect 
to the base functions, where the coefficients depend only on the internal variables. Since 
the base functions are polynomials, we are able to derive easily the generalized radial 
equations satisfied by the coefficients explicitly in Sec. IV. The permutation symmetry 
for the total wave function when some or all of the particles in the system are identical 
particles is discussed in Sec. V. In Sec. VI we will derive the generalized radial equations 
in a general case where the Jacobi coordinate vectors Rj (orthogonal vectors) are replaced 
by arbitrary coordinate vectors r^j in the center-of-mass frame (nonorthogonal vectors). 
In Sec. VII we will discuss a physical application of our approach. Some conclusions are 
given in Sec. VIII. 

II. SEPARATION OF MOTION OF CENTER-OF-MASS 

For a quantum A^-body system, we denote the position vectors and the masses of 
particles by and by m^, k = 1,2, . . . , N, respectively. M = J2k is the total mass. 
The Schrodinger equation for the A^-body system is 
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-- Y m,' Ar, ^ + V^ = i?vl/, (6) 

^ k=i 

where Ar^, is the Laplace operator with respect to the position vector r^, and \^ is a pair 
potential, depending upon the distance of each pair of particles, \rj — r^l- 

Now, we replace the position vectors by the Jacobi coordinate vectors Rj: 

Ro = M-'/^ Y mkVk, R, = (^^^) 

N 

l<J<(iV-l), M, = Ymk, M, = M, (7) 

k=j 

where Rq describes the position of the center of mass, Ri describes the mass-weighted 
separation from the first particle to the center of mass of the remaining particles, R2 
describes the mass-weighted separation from the second particle to the center of mass 
of the remaining N — 2 particles, and so on. The mass-weighted factors in front of the 
formulas for Rj are determined by the condition 

N N-1 
k=l j=0 

where an additional factor \/M is included in Rj for convenience. One may determine 
the factors one by one from the following schemes. In the center-of-mass frame, if the 
first j — 1 particles are located at the origin and the last N ~ j particles coincide with 
each other, the factor in front of R^ is determined by 

N 

Tj+l = = ■ ■ ■ = Fiv = -mjTj/Mj+i, Yl ^k^l = R^. (8) 

k=j 

A straightforward calculation by replacement of variables shows that the Laplace op- 
erator in Eq. (6) and the orbital angular momentum operator L are directly expressed in 
Rj : 

AT N-1 

k=l j=0 

N N-1 

L = -ihY rfc X Vrfe = -ifi Y Ri ^ Vr,, (9) 

k=l j=0 

In the center-of-mass frame, Rq = 0. Since the Laplace operator does not contain 
mixed derivative terms, the Jacobi coordinate vectors are also called the orthogonal vec- 
tors The Laplace operator obviously has the symmetry of the 0(3A^ — 3) group with 



respect to (3A^ — 3) components of (A^ — 1) Jacobi coordinate vectors. The 0(3A^ — 3) 
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group contains a subgroup SO(3)xO(A^ — 1), where S0(3) is the usual rotation group. 
The space inversion and the different definitions for the Jacobi coordinate vectors in the 
so-called Jacobi tree can be obtained by 0(A^ — 1) transformations. For the system 
of identical particles, the permutation group among particles is also a subgroup of the 
0(A^ — 1) group. As a matter of fact, after the transposition [k, k + 1) between the kth 
and the {k + l)th particles, the new Jacobi coordinate vectors, denoted by R^-, satisfy 

when j k or k + 1, 

1/2 r AT 



Mk (mfc + Mk+2) 



j=k+2 



Rfc cos6'fc + Rfc+i sin Ok, 

1/2 



{nik + Mk+2) Mk+2 

= Rfc sin 6k + Rfc+i cos 6k, 



N 



Mk+2rk - J2 

j=k+2 



cos6'i; 



mkUik+i 



Mk+i {ruk + Mk+2) 



1/2 



sin 6k 



MkMk+2 



Mk+i {rUk + Mk+2) 



1/2 



(10) 



This is obviously an 0(A^ — 1) transformation. For a system of identical particles, cos 6k 
{N-k)-\ 

It is easy to obtain the inverse transformation of Eq. (7): 



.1/2 



mjMj 



R, 



rrijMj^ 



1/2 



Ri 



-E 

k=l 

E 

i=k+l 



rrik 



MkMk+i 



MiM,+i 



1/2 



Rfc + M'"^Yio, 



- 1/2 


Mk 


Ri — 




_mfcMfc+i_ 



.1/2 



Rfc 



(11) 

(12) 



Thus, the potential is a function of Rj ■ Rfc. 

The Jacobi coordinate vectors Rj are invariant in translation and constitute a complete 
set of the coordinate vectors in the center-of-mass frame. If a complete set of arbitrary 
coordinate vectors Vcj is chosen to replace the Jacobi coordinate vectors. 



^cj ^ ' Rfc-Dfcj, 
fc=i 



det D^O, 



(13) 



where Dkj are functions of the masses rrij, the Laplace operator and the angular momen- 
tum operator become 

N-l N-1 

Sjk = DtjDtk, L = -i^ ^ r^j X Vr,,- (14) 
t=i j=i 

A typical example is v^j = rj — M~^/^Ro [see Eq. (11)]. When the matrix S is not diagonal, 
Tcj are called the nonorthogonal vectors |I5|JT6[] . We will not discuss the nonorthogonal 
vectors until Sec. VI. 



Af-l 

A = E ^.^v 

j,k=l 
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III. BASE FUNCTIONS WITH THE GIVEN ANGULAR MOMENTUM 



Because of the spherical symmetry, the angular momentum is conserved. We are going 
to discuss the wave functions with the given angular momentum and parity. From the 
given form (9), the eigenfunctions of the angular momentum are the homogeneous 
polynomials in the components Rji, of the Jacobi coordinate vectors Rj. 

For a quantum two-body system, there is only one Jacobi coordinate vector R = r, and 
the eigenfunction of the angular momentum is the spherical harmonic function Y^{6,(p). 
What is the generalization of the spherical harmonic function for a quantum A^-body 
system? A naive idea for generalization is to introduce the Euler angles, as was done by 
Wigner [ll,^-^. Is it necessary to introduce angular variables in the eigenfunction of the 
angular momentum? 

As is well known, the harmonic polynomial 3^^(r) = r^Y^{6,(f) is a homogeneous 
polynomial of degree £ in the components of r, which satisfies the Laplace equation as 
well as the eigenequation of the angular momentum. It does not contain angular variables 
explicitly. The number of linearly independent homogeneous polynomials of degree i in 
the components of r is 

t 

= Y.i^-s + l) = {i+l){i + 2)/2. 

s=0 

The number of homogeneous polynomials that can be expressed as a product of r ■ r and a 
homogeneous polynomial of degree {i—2) is N{i—2). Because N{i)—N{i—2) = 2i+l, the 
remaining homogeneous polynomials of degree i are nothing but the harmonic polynomials 

For a quantum three-body system there are two Jacobi coordinate vectors Ri and R2 
and three internal variables in the center-of-mass frame: 

ei = Ri-Ri, 6 = Ri-R2, r/2 = R2-R2. (15) 

The internal variables are invariant in the global rotation and the space inversion of the 
system. We are going to construct base functions for angular momentum that do not con- 
tain a function of the internal variables as a multiplying factor, because the factor should 
be incorporated into the generalized radial functions. The number of linearly independent 
homogeneous polynomials of degree i in the components of the Jacobi coordinate vectors 

is M{ey. 

M{i) = ^{i +!)(£ + 2){i + 3)(£ + 4)(£ + 5). 
5! 

The number of the homogeneous polynomials of degree i that do not contain a function 
of internal variables ClS db factor is 

K{i) = M{i) - 3M{i - 2) + 3M{i - 4) - M{i - 6) = Af + 2, i> 1. 
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On the other hand, the wave function with angular momentum i can be obtained from 
y^(Ri)3^^''(R2) by use of the Clebsch-Gordan coefficients {q^m^i — q.^i — m\L^[i) 

3^i^,(Ri, R2) = E yi{^^)ylZli^^) {q, m, £ - g, /i - m|L, /i). (16) 

m 

3^^^(Ri,R2) is a homogeneous polynomial of degree £ in the components of the Jacobi 
coordinate vectors Rj. Simultaneously, it is the common eigenfunction of L^, and the 
space inversion with eigenvalues L{L + 1), and (—1)^, respectively. When fi = L = i 
(0 < g < £) and (£ - 1) (1 < g < £ - 1), we have 

V^'^m R^-r 1^J [(2g + l)!(2^-2g + l)!]^/M rP ^.ffi- ^^"'^ 

(Ri, R2j - (-1) < g!(£ - g)!2^+% I ^ ^^'^ ^ ^^'^ 

^ { (2i±||^}" {(, - 1),, - , ^ 

It is evident that these expressions do not contain a function of the internal variables as 
a factor, neither do their partners with smaller /i due to the spherical symmetry. The 
number of those eigenfunctions is 

(2£ +!)(£+!) + (2£ - 1)(£ - 1) = + 2 = K{1), i>\. 

That is, any of the remaining eigenfunctions 3^^'^(Ri, R2) with L < l — \ can be expressed 
as a combination, where each term is a product of a function of the internal variables and 
a homogeneous polynomial of degree less than i For example, 

yll{^^, R2) = -^6, J^o'o (Ri, R2) = 1^ {3^2' - ^1^2} , 



J^s's (Rl, R2) = {r72(i?lx + iRxyf + ei(^2x + lR2yf 

In other words, any eigenfunction with angular momentum ^ is a combination of those 
homogeneous polynomials 3^£^(Ri,R2) and 3^£^'''^''^(Ri, R2), where the combinative co- 
efficients are functions of the internal variables. Since the normalization factor can be 
ignored, we rewrite 3^j^^+^^''(Ri, R2) in a simplified form as Qg'^(Ri,R2) by removing a 
constant factor 

Qf(Ri,R2)= A<g<£, A = 0,1 

X = Rix^iR\y, Y = R2x + iR2y, Z = X — RizY- (18) 
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Note that 

gf (Ri, R2) = Qi'll^\Ri, R2)Z. (19) 

(5g^(Ri,R2); called the generalized harmonic polynomial, is a homogeneous polynomial 
of degree {£ + A) in the components of the Jacobi coordinate vectors. It is the common 
eigenfunction of L^, L^, h^^^, Ar^, Arj, VRi ■ VR2) ^^e space inversion with 
the eigenvalues £{£ + 1), £, q{q + 1), {i - q + X){i - g + A + 1), 0, 0, 0, and (-1)^+^, 
respectively, where L^^ (Lrj) is the square of the partial angular momentum, and Ar^ 
(ArJ is the Laplace operator with respect to the Jacobi coordinate vector Ri (R2) [see 
Eq. (9)]. Any wave function with the given angular momentum i and the parity (— 1)^+^ 
can be expressed as follows: 

vl/f (Ri,R2) = Yl V^f (6,6,r/2)Qf (Ri,R2), A = 0,1. (20) 

That is, for a three-body system the generalized harmonic polynomials Qg^(R.i, R2) con- 
stitute a complete set of base functions with angular momentum i and parity (— 1)^+'*'. 
Only £ -|- 1 — A partial angular momentum states are involved in constructing a function 
with angular momentum i and parity (—1)^+^, and the contributions from the infinite 
number of remaining partial angular momentum states are incorporated into those of the 
radial functions. Substituting Eq. (20) into the Schrodinger equations (6) and (9), one is 
able to easily derive the generalized radial equations for the generalized radial functions 



^ nei,6,^2) Mm 



AV^f + 4g%V^f + 4{i-q + A)9,,^f + 2{q - A)%V^J^ + 2{i - g)%^J^i 

2 



A^^ei, 6, V2) = {46^1, + ^V2d',, + 6 (% + d,,) 

+ (6 + m) dl + 46 (% + % } (^1, 6, V2), 

X<q<i, A = 0,1. (21) 

For a quantum A^-body system, there are (A^ — 1) Jacobi coordinate vectors. We 
arbitrarily choose two Jacobi coordinate vectors, say Ri and R2. We fix the body-fixed 
frame such that Ri is parallel with its Z-axis, and R2 is located in its XZ plane with a 
non- negative x-component. Define (3A^ — 6) internal variables, which are invariant in the 
global rotation of the system: 

Cj = Rj ■ Rl; Vj ~ Rj ' -^-2; — Rj ' (-R-l ^ 1^2) ; 

1<J<(A^-1), m = ^2, Ci = C2 = 0. (22) 
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It is worth mentioning that and rjj have even parity, but Q has odd parity. From them 
we have 

Qj = (Ri X Kj) ■ (Ri X R2) = CiT]j - ^2^j, 
ujj = (R2 xRj) ■ (Ri X R2) = ^2Vj - V2^j, 

Q^=uj2 = 0, Q2 = -uJi = (Ri X Ra)^ . (23) 
Due to our choice of the body-fixed frame, the components of Ri and R2 in the frame 



are 



(0,0, 



/2 



and 



(^2/6 



a/2 



0,6^1 



■1/2 



respectively. From Eq. (22) we are able to 
express all the components R^;, of the Jacobi coordinate vectors Hj in the body-fixed 
frame by the internal variables: 



nj (6f^2) 



A/2 



R 



QVI2 



'1/2 



R 



-1/2 



(24) 



The formulas (24) also hold for j = 1 and 2. The volume element of the configuration 
space can be calculated from the Jacobi determinant by replacement of variables: 



N-l 



1 



Y[ dRjxdRjydRjz = -^l sin I3dad(3d'yd^id^2dri2 Y\_ d^jdrjjdQ. 



(25) 



The ranges of definition of the Euler angles are well known, the ranges of definition of ^1 
and r]2 are (0, 00) and the ranges of definition of the remaining variables are (—00, 00). 
Furthermore, 

Rj ■ Rfc = ^2^ i^jVk - ^j^k + CjCfc) • (26) 

It is easy to see from Eqs. (12) and (26) that the potential is a function of only the 
internal variables. Since Ri and R2 determine the body-fixed frame completely, it also 
can be seen from Eq. (24) that each of the components of the Jacobi coordinate vectors Rj 
can be expressed as a linear combination of R16 and Rah with the coefficients depending 
on the internal variables. In fact, denote the rotation transforming the center-of-mass 
frame to the body-fixed frame by -R(a,/3, 7) with three Euler angles 



/2(a,/3,7) 



f5 Q( Cry I S ry ^ Oi ^ j3 ^ '^f l ^ct ^ (X^ (3 

\ ~SpC^ SpS^ Cf3 J 



(27) 



where = cos a, 
(24) and (27) that 



sin a, and so on. It is straightforward to obtain from Eqs. (18) 



Y = R2, + iR 
a/2^ 



X — Rlx + iRiy 

.1/2 



2y 



Rlz = Cl^Cp, R2z = - {^2/^1)^^^ Sf^Cy + ^2^1 



1/2 



C/3, 



Z = {Rlx + iRiy) R2Z — Rlz {.R2X + iR 



2y) 



{2i 
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Rjx + iRjy = ^2 ^ + ~ ^CjZ} , 

{Rjx + iRjy) Rkz — Rjz {Rkx + iRky) = ^2 ""^ iVjCk — VkCj) ^ 

-I (Oa - 60) y + (O^fc - ^km) z} . (29) 

Therefore, each harmonic polynomial y^(Rj) can be expressed as a combination of 
Qq^(Ri,H2) with the coefficients depending on the internal variables. This means that 
the generalized harmonic polynomials (^^^(Ri, R2) given in Eq. (18) do constitute a com- 
plete set of independent base functions with the given angular momentum i for a quantum 
A^-body system, just like they do for a quantum three-body system. 

Because this conclusion plays a key role in separating the global rotational degrees of 
freedom from the internal ones in the quantum A^-body system, we are going to prove it 
by another method. From Eq. (28) we have 

e'^'sp = ^i^^^X, e*" (c^ + icfss^) = -Vl^^^'^Z, 

e'" (c^c-, + ts^) = -6 (6^2)"'/' X + {ii/n^f^ F, (30) 

= 7]2X^ - 2^2XY + iiY^. (31) 

That is, a homogeneous polynomial of degree I in three variables e^'^sp, e*" {cpc^ + is^), 
and e*" (c^ + icps^) can be expanded with respect to (^^'^(Ri, R2) where the coefficients 
only depend on the internal variables 0, Vji Q. On the other hand, the Wigner 
/^function is 



e*" (c^ + ^C/3S^) ± e*" {cpc^ + is^) , 

(32) 

where m > 0. Therefore, D£^(a, /5, 7)* can be expanded with respect to Qg'^(Ri,R2), 
where the coefficients depend only on the internal variables. Z)^.^(a, /?, 7)* constitute 
a complete set of independent base functions with the angular momentum £, so do 
Qg'^(Ri,R2). Now, we come to the theorem. 

Theorem. Any function \E'£'^(Ri, . . . , R^v-i) with the angular momentum £ and the parity 
(— 1)^+'^ in a quantum A^-body system can be expanded with respect to the generalized 
harmonic polynomials (^^^(Ri, R2) with the coefficients rj, () depending on (3A^ — 6) 

internal variables: 

vi>f (Ri, . . . , R^_o = E E i'lrU^ V, OQf (Ri, R2), 

r=0 q=-r 



{2i)\ 



{e + m)\{e-m) 



1/Z 



bp 



(1 ± C/3) 



-1^ 



(2^)! 
_(£ + m)!(£-m) 
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i^fM^ V, = i^lri^l, ^N-l, V2,---, VN-1, Cs, • • • , Cn-i), 



C(e,r/,-C) = (-l)'->^^(e,r/,C), (33) 



where the last equahty means that the parity of '?/'q:^(^, r], () is (—1)'^ ■ 



IV. THE GENERALIZED RADIAL EQUATIONS 

From the theorem above, the set of Qg^(Tli, R2), just hke the set of Df^^^a, P, 7)*, is a 
complete set of base functions with angular momentum i in the quantum A^-body system. 
Each function with the angular momentum i in the system can be expanded like Eq. (33) 
or Eq. (2). However, Eq. (33) has two important characteristics, which make it easier to 
derive the generalized radial equations. One is that the generalized harmonic polynomial 
Qg'^(Ri,R2) is a homogeneous polynomial in the components of two Jacobi coordinate 
vectors Ri and R2, where the Euler angles do not appear explicitly. The other is the well 
chosen internal variables (22), where the internal variables Q have odd parity. It is due to 
the existence of Q that (5g°(Ri, R2) and (^^^(Ri, R2) appear together in the expansion of 
the wave function. By comparison, all internal variables in a quantum three-body system 
have even parity {(j = 0) so that in the expansion (20) of a wave function with a given 
parity only the base functions with the same parity appear 0J^,^,2S . 



Because of these two characteristics, it is easy to derive the generalized radial equations 
by substituting Eq. (33) into the Schrodinger equation (6) with the Laplace operator (9). 
The main calculation in the derivation is to apply the Laplace operator (9) to the function 
\l/£'^(Ri, . . . ,RAr_i) in Eq. (33). The calculation consists of three parts. The first is to 
apply the Laplace operator to the generalized radial functions iplri^yVX)'- 

A^J^(e, V. = {4ei5|, + ^V2dl + (6 + V2) dl + 46 (% + d^, 

N-1 

+ 6 (% +d,,)+Y. [^idl + V2dl + n,dl + 269^,9,^. 

+ 4 [iA. + CA) % + 4 {nA, + CA) K + 2 (r^A, + iA^ 
+ ^2^ [(^j^k - ujj^k + oCfc) {d(Ak + dvAk) - 2 (^iCfe - ^fcCi) hAk 

j,k=3 

+ 2 (^-a - ^kQ) dr,A. + i^j^k + ujjujk + UjCk + V2QCk) 9(A.] } V^Jr (e, V, 0- (34) 
The second is to apply it to the generalized harmonic polynomials Qi'^{T{.i, R2). This part 



'9 



is vanishing because (Ri,R2) satisfies the Laplace equation. The last is the mixed 
application 

N-1 



2 \ 2Ri + R2 + X: [{d^^i^lr) R. + (^C.V^J') (R2 X R,)] \ ■ VR.a 
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N-l 



In terms of Eqs. (18) and (29) we obtain: 

Ri • VRiQf = qQl\ R2 • VR.Qf ^(e-q + r) gf , 

R2 • VRiQf ={i-q + l) gf_i, Ri • VR.Qf = {q-T + l) Qf+i, 

R,- • Vr,< = ' {-^.-^Qf + Ji,- (^ - ? + 1) gf_i - iOQf } , 
R. ■ Vr.< = ^2^' {-^.(g + Will + - q) Qf - <jQl\i} , 

R,- • VriQ? = ^^2"' {-mCjq'Qf + ^UA^q - 1)(^ - ? + l^-i 
- i^iCA^ -q + 2){e-q + l)Qf_, - UjqQf + n^(e-q + l)Qf_^} , 

Ri • Vn.Qf = {-mCAq + ^)qQf+i + ^60^(2^ - 2g + i)gf 
- - q + - ^.gQf+i + - q + } , 
(R2 X R,) ■ VR.Qf = ^^2"' {^20gQf - 60 (^ - g + 1) Qf-i - } , 

(R,- X Ri) • VR.Qf = ^^2-' {-60(g + 1)Q;% + CiC, (^ - q) Qf + ^^JQl\l} , 
(R2 X R,) • VRiQf = ^^2 ' {-ir/2a;,g2gf + ^6c^,(2g - 1) (^ - g + 1) Qf_i 
- z6^,(£ -q + 2){£-q + l)Qf_, + rj^QqQl' - 60 (i-q + l) Qf-i} , 
{Rj X Ri) • VR.Qf = ^^2 ' {iV2^j{q + l)qQf+i - ^6^^,? ^ -2q + l) Q[ 

+ zei^(^ - g + i)'gf-i - 60gQf+i + 60 (^ - g + 1) } • 

Now, the generalized radial equations are 



A< + 4 {g^a + (£ - g)a,J + 2qd^,^l^_,), + 2(£ - g)%4\i)o 



//A 



Af-1 



+ E 2f^2-' { 



iA 



J=3 



+ ^■(^ - q)dni + mCjqdci + 60 (^ " g)^^] ^go 

- g [ujd^. + 60^0] ^S-i)o + - [^Aj - 605cJ ^S+i)o 
- iwiq - 1) [O^r,^ - ^A^ ^{q-l)l 



-iq 



^20?%- - 60(2^ - 2g + 1)9,, + r72a;,-g9c, + 6^^i(2^ - 2g + l)ac,J 

+ ^{i - q) [60(2g + 1)% - 60(^ - qA. + 6^i(2g + i)ac, + 6%(^ - ?)5g] V^g+m 

-m -q){^-q-l) [0% + ^i^c] C+2)i} = - (2/^') - V] i^ll (36a) 



(35) 



A^l>l + 4 {g% + (£ - g + 1)9,,} .^^^ + 2(g - lA^Un + " 
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Af-1 



j=3 



{q - 1) [u.d,^ + 605c.] V^g-Di + - [^.5«. - 605g] 



- ^^j^o 



V'S-Do - i [Cjd^^ + ^A,] V'S} = - (2/^') - V] V'S, (366) 



where A'^^:^ was given in Eq. (34). When A'"=3, Eq. (36) reduces to Eq. (21), where, 
because all internal variables have even parity, the generalized radial functions ipq'^{C, ?7, C) 
with T have to be vanishing. 



V. PERMUTATION PROPERTY OF WAVE FUNCTIONS 

When some or all particles in a quantum A'"-body system are identical particles, one 
has to consider the permutation property of the spatial wave function, which depends on 
the total spin of identical particles. Since the spatial wave function ^|'^(Ri, . . . , Riv-i) 
is expanded with respect to the base functions (5g^(Ri,R2), we need to study only the 
property of (5^^(Ri, R2) in the transposition {k, A; + 1) between two neighboring particles. 
The transformation property of the Jacobi coordinate vectors Rj in the transposition 
{k,k + 1) was given in Eq. (10). Therefore, the base function (5^^(Ri,R2) remains 
invariant in the transposition (A;, k + 1) with A; > 3. In the following we are going to study 
the transformation property of (^^^(Ri, R2) in the transpositions (1, 2) and (2, 3). Denote 
by Pi and P2 the transformation operators for the base function in the transpositions 
(1,2) and (2,3), respectively. In the following formulas wc neglect the argument Ri,R2 
in Q^^(Ri,R2) and briefly denote sin 6*1, sin 6*2, cos 6*1, and cos 6*2 by si, S2, ci, and C2 for 
simplicity. 

A. Transposition (1,2) 



where 

n {q ~ n)\{l - p - n)\n\{n + p - q)\' ^ ' 

Because of Eq. (19) and PiZ = —Z, we obtain 

p.Qf = -i: (39) 

p=l 

B. Transposition (2, 3) 
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E QfDlf{92)+^ E Q?^Jf(^2), (40) 



P=? p=q+l 



where 



r>^(2) _ (g + n)!(£-g-n)!4 ' ^ ( 1 /_1^n+m^2m 



rt 



— q — 2m \ ( 'm \ ( m — t 
n-r ]\t]\r-2t 



n^(3) _ (g + n -!)!(£- g - nyJf ^ / i-q \ , .n+m^2m+i 



xE 



- q — 2m — 1\ I m \ I m — t 
n-r -I /\t/\r-2t 



X (26)'--'* ujt'-' (^3 - 1^202/^2)'-^-"-'"^+'' , (41) 

where the combinatorics ^ 5 ^ ~ ^[^ — ranges of the summation indices m, 
r and i are determined by the conditions that those combinatorics are not vanishing. 

= i«2C3^^2-^ E - Q?-i(^ - P + 1)6} (^2) 

p=q 

+ i {s^n, - 02^2) ^^2"' E {Q>ip - 1)^72 - 2Qf_,{p -i)ii-p + 1)6 

p=q+l 

+Q'U£ -P + 2){e-p + 1)6} ^(tliS-i) (^2) 
+ (.2^13 - C2Q2) E QpK-V-i)(^^) 

p=q 

-«2C3^^2-' E {<(p-i)6-Q^i(^-P + i)6}/^(t;)(ti)(^2). (42) 

p=9+l 

In real calculations the cases with the small angular momentum may be more inter- 
esting. In the following we explicitly list the above formulas for £ — 1 and 2 (the formulas 
for the case with £ — are trivial) : 
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PiQi 



PiQf 



20 



-2ciSiQ 



20 



ci-si]Qf + 2c,s^Q-f, 



20 



^21 



21 



->21 



•>21 



ill 



P^Qf = Qf, P^Qf = 82^-^ [-^OJsQf + (Q3 - ^^2C2/S2) Qf - iCsQ'. 

in^' { (a;| - mC!) Qf + [-UJ3 (1^3 - f^2C2/s2) + 6C|] QT 



P2Q, 

+ 



20 




(Qa - Q2C2/S2)' - ClCi] Ql' + i^sCsQf - <3 (^^3 - ^^2C2/S2) Qf} , 

P2Qf = 82^2' {^26(3^2° - ^CiCsQf + (^^3 - ^202/52) Qf} , 

P2Qf = Slnf {-2lCs [^2003 + m {^3 - ^^2C2/S2)] Qf 
+ <3 [6^3 + 36 (^^3 - ^^2C2/S2)] QT - 4i6C3 (^^3 - ^^202/52) 

+ [6C| - ^3 (^3 - ^^2C2/S2)] Ql' + [(^3 - ^^2C2/52)' - 6C|] QI'} ■ 



VI. NONORTHOGONAL VECTORS 



Now, we turn to the general case where arbitrary coordinate vectors in the center-of- 
mass frame [see Eq. (13)] are used to replace the Jacobi coordinate vectors Rj. In this case 
the Laplace operator contains mixed derivative terms [see Eq. (14)]. All the conclusions 
in Sec. Ill hold for the present case except that the Jacobi coordinate vectors should 
be replaced with the coordinate vectors r^j and the volume element of the configuration 
space (25) changes due to the linear transformation (13). In particular, the generalized 
harmonic polynomial (5g^(Ri,R2) now becomes (^^^(rci, rc2), where and rc2 are two 
arbitrarily chosen coordinate vectors. 

Any function ^|'^(rci, . . . , rc(jv-i)) with angular momentum i and parity (—1)^+^ in a 
quantum A^-body system can be expanded with respect to (5g^(rci, i"c2) with the coeffi- 
cients ■05t(C' V: C) depending on (3A^ — 6) invariant variables 

*f (ra, . . . , r,(^_i)) = ^ ^ Vj^(e, 77, OQf (rci, r,2), (43) 

T=0 q=T 

where, instead of Eq. (22), the internal variables ^j, r)j, and Q are redefined as 
~ ^cj ■ Tel, Tjj = Tcj ■ rc2, Cj = ^cj ■ (rci X rc2) , 
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i<i<(A^-i), 771 = 6, Ci = C2 



0. 



(44) 



As in the case with the orthogonal vectors, the main calculation in deriving the gen- 
eralized radial equations in the present case is to apply the Laplace operator (14) to 
the function \l>^'^(rci, . . . , rc(Ar_i)) in Eq. (43). Similarly, the calculation consists of three 
parts, and the second part is vanishing. But, the first part [see Eq. (34)] becomes 

AVj^(e, V, C) = {(4^116^1, + 4^22772^4) + {Sum + ^226 + 2^126) dl 

+4 (5ii6 + Su^i) + 4 (^226 + ^12^2) 



N-1 



+ E {4 (Su^j + -51,^1) + 4 (^227;, + S^.m) d,,d,^ + 4 {Su^j + 5i,6) d^A 
3=3 

+ 2 {S22i, + SuVj + <5i,7;2 + 52,6) + 4C, (-Sn9^, + ^22^,, + ^12%) d^. 

N-l 

+ 6 + S2A,) } r/, C) + E {(2-51.^. + ^,.6) d^A, 

j,k=3 

+ i'^S2jVk + Sjkm) dnjdr,^ + 2 (SikVj + S2j^k + -^^^6) d^^dr,^ + 2 {SijCk + SikCj) ^^.-^a 
+ 2 {S2jCk + ^2^0) 5,,% + i2S^,uJk - 252,i], + ^.-.fis) d^^^, } ^' C) 



N-l 



+ E ^2 ^ {{^jVk - ^^jCk + CjCk) [Sud^.d^^ + S22dr,Ar,k + '^S^d^^, 
j,k=3 



Vk 



+ 2 [Sii (ujCk - uJkCj) + S12 {^jCk - ^kCj)] d^jdi;^ 

+ 2 [5'22 {^jCk - ^kCj) - S12 {ijjjCk - ^fcCi)] ^Ty.^Cfe + ['S'li i.^j'^k + mCjCk) 

+ ^22 {^j^k + Cioa) - 2-5i2 {uj^k + 6oa)] V'Jrie, r;, 0- 

The last part contains the mixed application 



(45) 



N-l N-l 

2 E E s,k vr., < • vr.,g: 



9 



2^r 



cl 



25n% + 5i2% + E -5^% 

j=3 



rc2 



Af-l 



3=^ 



N-l 



+ E \^<^3 (-^ll^^i + -51257,^) + (rc2 X Tcj) ^ii^^. + {Vcj X Ta) 5i25cj 



j=3 



+ (rei X r,2) ^y^c,] } • Vr.,gf + 2 



Af-l 



2^12% + -522% + E '^2,'95^ 
j=3 
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+ rc2 



N-l 



j=3 



N~l 



+ (rc2 X Tcj) Si2d^^ + (r,j x r,i) ^22% + (rci x rc2) 52^9^^ } i^l^ ■ Vrc2<3f • 



(46) 



In addition to the formulas (35), where Rj should be replaced with Tcj, we also need the 
following formulas 



(re2 X r,,) • Vr..gf = {V2CA1 + - 60(^ - g)Qf - , 

(rc2 X re,) • Vre.gf = ^^2"' {-^^i^ji^ ' Q + l)'Q?-i + ^^2UJjq i2i -2q + l) gf 
- irj^cojiq + l)?gj% - - q + l)Qf + V2CjqQl+i} , 

(re, X rei) ■ Vr.,gf = ^^2"' - q + - Q + 2)Qf_2 

- (^ - g + 1) (2g - l)gf_i + ^ry2^•g'gf + UA^ - q + l)Qf-i - U.qQf} , 



(rcl X Tea) • VrciQf = -iQ^q, 
(rel X re2) ■ Vr^^gf = -iQ^+l, 

(rei X re2)-Vr.igf = -^772? + ^6 " 1) (2? - l)gf_i - - ? + 2) (^ - ? + l)gf_2 , 
(rel X re2) • Vr..gf = ^6^(2^ " 2? + l)gf - - g + l)'gf-i - iV2q(q + 1)Q% (47) 
Finally, we obtain the generalized radial equations as follows 



ex 



AV'S + 2 [2Snqd^, + iS^d^, + 2(£ - q)S22dr„] V'.o 
+2q {2Sud^, + ^22%) V'g-Do + 2(£ - g) (-Sn^^, + 23^0,,) V^g+i)o 

+ ^ 21^2 1 { [(£ - q) (n^S^2d^^ + 1^,522^,^. + n2S,2d,^ + 60^22^0 



j=3 



//A 



+q {^l2Sji - ujjSii) df:^ - ujjqSudr,^ + Cj {ri2qSu - ^2^812) <9^^-J ip^Q 

+q [{^2Sj2 - UJjSu) d^. - UJjS22dr,j + Q {^Su " 6'S'22) ^Ci] -0(^-1)0 

+{£ - q) [Sn^A, + iSji^2 + Si2nj) 9,. + (^1-512 - 6-5ii) ac,] V(,+i)o 

-zrj2q{q - 1) {Q (-Si2% + -522^,^ + (cu,,Si2 - ^-^22 + ^^25',-2) d^^ } i^ll,^, 
-iq r)2q {CjSiAj + CjSudrjj + oUjSndQ - ^jSi2dQ + 9.2SjidQ^ 
-6(2£ - 2g + 1) {QSi2d^^ + QS22d,^ + ^28^20^, - ^,8228^^ + o^,^i25o)] < 
+i(£ - q) [{2q + 1)^2 (O^ii^^, + 0'5i2\. + ^^2S',iac, + UjS^A. - iljSudc,) 
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'/A 



<l{^ - g)(^ - g - 1) [O (^11% + 5i29,J + (cU,-Sn - ^^,512 + ^^25',l) 

= -(2A2)[i5;-V]V^3, (48a) 
A7/>2 + 2 {2g5n% + + l)5i2% + 2(£ - g + 1)5229,J 
+2(g - 1) (2^12% + ^22%) + 2(£ - g) (5n% + 2S,2d,,) 



iV-l 



J=3 



-'(9-1)0 



//A 



+ [g (il25'ji - LOjSii) d^^ - ujjqSi2dr,j + Q {r]2qSu - ^2(^ + l)5'i2) 
+ {i-q + 1) {n.Sud^^ + ^28,28.,^ + n^S22d^^ + U,S22dc,)] t 



/A 
9I 



+(g 



) [(^25'j2 - UJjSu) - UJjS22dr,^ + (j {V2S12 - ^2822) ^C.J V^^-l)! 

+ (£ - g) [njSnd^^ + {^28^, + n,Su) + Q (-6^11 + 6^12) 5^,] 



VII. PHYSICAL APPLICATION 

In a quantum A^-body system, any function with angular momentum i and par- 
ity (—1)^^'^ can be expanded with respect to the generahzed harmonic polynomials 
(5g^(Ri,R2)5 where the coefficients, called the generalized radial functions, depend only 
on (3A^ — 6) internal variables. Since Qg'^(Ri,R2) is a homogeneous polynomial in the 
components of the Jacobi coordinate vectors and a solution of the Laplace equation, we 
have derived the generalized radial equations easily. That the rotational variables (the 
Euler angles) do not involve in either the generalized radial functions or the equations will 
greatly decrease the amount of calculation in solving the Schrodinger equation numeri- 
cally for the A^-body system. As a ffist step, we applied this approach to the calculation 
of the energy levels of a helium atom and a positronium ion |p9H3T[| . In the following we 
sketch the method and give some more calculation results. 

Once the generalized radial equations have been derived, one may choose any other 
complete set of internal variables to simplify the calculation. The generalized radial 
equations for the new variables can easily be obtained by replacement of variables. In 
a Coulombic three-body system, such as a helium atom, we choose the hyperradius p 



20 



and two dimensionless rj and ( as the internal variables, so as to make the potential a 
meromorphic function: 



P=(r+r) , r, = ^ ^, C = ^ - + — 49 

where ri denotes the position vector of the helium nucleus and r2 and the position 
vectors of two electrons. After expanding the wave function as a Taylor series with 
respect to rj and ( up to the order n, we obtain an ordinary differential matrix equation 
for the coefficients R{p). In the real calculation, we calculate the propagating matrix F{p) 
and its inverse matrix G{p) by the Taylor series method instead of the function R{p) in 
order to avoid the logarithmic singularities at p = in the forms of p°(lnp)^ |]32|-p4|: 

= F{p)Rip). (50) 

We are able to obtain the energy levels of a helium atom in the different spectra '^^+^L'^^"^ 
with high accuracy by choosing n = 10 (for the S and P states) or n = 9 (for the D 
states) due to the fast convergence of the series, where S is the total spin of two electrons, 
e(o) describes the parity, and L = S, P and D for the angular momentum states. In 
order to compare our calculation results with those by the variational methods where the 
nucleus mass is usually assumed to be infinite, we also calculate the energy level with a 
large mass ratio M of the nucleus to the electron (M = 10^°). Both calculated results 
are listed in Table 1 for comparison. Much fewer terms in the truncated series are taken 



in our calculation than those in the hyperspherical harmonic function method and in 
the variational methods . 

Table 1 



VIII. CONCLUSIONS 



For a quantum A^-body system we have found a complete set of independent base 
functions Q^^(Ri,R2) for the given angular momentum and parity. Any function with 
angular momentum i and parity (—1)^+'^ in the system can be expanded with respect to 
the {2i + 1) generalized harmonic polynomials Qq^(Ri, R2); where the combinative coeffi- 
cients are functions of the (3A^ — 6) internal variables. We have established the generalized 
radial equations explicitly. They are simultaneous partial differential equations in the in- 
ternal variables. The number of both the functions and the equations is {2i + 1) when 
A^ > 4, and it becomes {£ + 1) or £ when A^ = 3, depending on the parity. Only a finite 
number of partial angular momentum states are involved in constructing the generalized 
harmonic polynomials (5^^(Ri, R2)- That is, the contributions from the remaining partial 
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angular momentum states have been incorporated into those from the generahzed radial 
functions. We have generalized the formulas to the case with nonorthogonal vectors. 

When establishing the body-fixed frame we fix it with two arbitrarily chosen Jacobi 
coordinate vectors Ri and R2. Those two vectors may be replaced with any other two 
coordinate vectors according to the characteristics of the physical problem under study. 

The choice of the complete set of base functions is not unique. However, the right 
choices of both the base functions and the internal variables play a key role in establishing 
the generalized radial equations. Those two choices are the main progress of the present 
paper in comparison with the previous work of Wigner [^] and Eckart [^. Once the 
generalized radial equations have been derived, one may choose any other complete set of 
the internal variables to simplify the calculation. The generalized radial equations for the 
new variables can easily be obtained by replacement of variables, just as we did in Sec. 
VII for the three-body system. 

The two features in this method, that the numbers of both functions i^lri^^VjC) 
equations are finite, and they depend on only (3A^ — 6) internal variables, are important 
for calculating the energy levels and wave functions in a quantum A^-body system. In 
fact, in the numerical experiments for a quantum three-body system, we calculated the 
lowest energy levels of a helium atom in P states and in D states [^| with the total 
spin one and zero, and some energy levels of a positronium negative ion |^ by a series 
expansion. Because three rotational variables are removed, many fewer terms have to be 
taken to achieve the same precision of energy as in other methods to truncate the series 
of partial angular momentum states. As the number of particles in the system increases, 
we believe that removing three independent variables related to the global rotation will 
greatly decrease the amount of calculation. 
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Table 1 Numerical calculation for the energy levels of a helium atom in atomic units 



Spectral Term 

2S+l^e(o) 


Our Results 


Variational Calculation [Q 

M ^ oo 


M = 7296.28 


M = 102" 


l^e 


2.9033046 


2.903724377034116 


2.9037243770341195 


S^e 


2.1749303 


2.1752293777 


2.1752293782 


Ipe 


0.5801748 


0.5802464725 


0.5802465t 


3pe 


0.7105002 


0.7103965 


0.710499^ 


Ipo 


2.1235456 


2.1238430778 


2.1238430865 


3po 


2.1328807 


2.133164187 


2.133164191 


lp)0 


0.5637256 


0.5638004 




3p)0 


0.5592482 


0.5593283 




3p)e 


2.0553230 


2.0555871 


2.0556363 


lj~)e 


2.0553055 


2.0555693 


2.0556207 



^ The calculation in Ref. |^| 
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